Effect of Interatomic Separation on Entanglement Dynamics in a Two- Atom 

Two-Mode Model 



K. Sinha0 N. 1. Cummingsjt] and B. L. 

Joint Quantum Institute and Department of Physics 
University of Maryland, College Park, Maryland 20742-4111, USA 
(Dated: August 6, 2011) 

We analyze the time evolution of quantum entanglement in a model consisting of two two-level 
atoms interacting with a two-mode electromagnetic field for a variety of initial states and interatomic 
separations. We study two specific atomic separations which give rise to symmetric atom- field 
couplings. For general atomic distances we consider a subset of initial states analytically, and then 
treat the more general situation numerically. We examine a variety of qualitative features such as 
entanglement sudden death, dynamical generation, protection, and transfer between subsystems. 
Our analysis shows a stark contrast in features of entanglement between the two special coupling 
schemes often considered; further, these features are uncharacteristic of those arising for general 
distances, due to the high degree of symmetry present in the special cases. The variety of behaviors 
in these two-mode cases suggest the importance of considering atomic separation carefully for any 
model where two atoms interact with a common field. 



I. INTRODUCTION 

Quantum entanglement has been extensively studied, 
both due to its fundamental significance in quantum the- 
ory [1, 2J and its utility as a resource for quantum com- 
munication and quantum information processing [3, 5j. 
Atomic physics offers a domain with sufficient control of 
the system and isolation from noise that it has been pos- 
sible to perform precision experiments on quantum en- 
tanglement [Sll^, and, therefore, it is also a productive 
target for theoretical study of the issue as well. 

One of the simplest scenarios for theoretical studies 
of the dynamics of entanglement between atoms is that 
of two atoms which are isolated from one another and 
interact with different electromagnetic fields. Studying 
this sort of model led to the discovery of the entangle- 
ment sudden death (SD) phenomenon [8HT2]. in which 
entanglement decays to zero in a finite time rather than 
asymptotically. SD has garnered significant interest be- 
cause it is both unintuitive and potentially undesirable. 

An alternative simple model for atom-field interaction 
in which to study entanglement dynamics is the model 
studied by Dicke [13j and Tavis and Cummings [14J, where 
one assumes all atoms are grouped in a sample whose 
size is small compared to the resonant wavelength (re- 
sulting in identical coupling to every atom). However, 
such a simplistic model may miss the variety of behav- 
ior that can result when the atoms are not confined to 
such a small sample. Entanglement dynamics have been 
shown to have a significant distance dependence when 
two atoms are interacting with a common field [15j. For 
atoms weakly interacting with a continuum of field modes 
in the Born-Markov approximation, it has been shown 
[16 - 19j that changing the atomic separation of two atoms 
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can affect whether there is SD and whether there is re- 
vival of entanglement, as well as modify the dynamical 
generation of entanglement; in short, the qualitative fea- 
tures are sensitive to the atomic spacing. At short inter- 
atomic distance non-Markovian effects associated with 
induced interactions between the atoms due to the quan- 
tum field become more pronounced whose qualitative be- 
havior varies greatly with different classes of initial states, 
as detailed in the study of [20j . So it is clear that distance 
dependence gives rise to a significant variety of behaviors, 
and some seemingly innocuous approximations can qual- 
itatively alter the entanglement dynamics in unintuitive 
ways pTl[22] . 

In view of this, we seek to study entanglement dynam- 
ics in a model that is sufficiently complex to manifest 
some of this variety of behavior yet simple enough that 
the dynamics may be understood in considerable detail 
and obtained with fewer approximations. The simplest 
model one may pursue along this line would be two two- 
level atoms coupled to a single mode with couplings that 
reflect the distance dependence; however, as discussed in 
Sec. [TTj one needs to include at least two field modes be- 
fore any non-trivial distance dependence can arise in the 
problem. So we will study such a model with two field 
modes here. We choose to have traveling wave modes so 
as to eliminate the effect of atomic center of mass position 
and focus on the relative separation between the atoms, 
this can be experimentally realized in the situation where 
one has two counter-propagating modes in a toroidal res- 
onator cavity ^38] . The motivation for working with 
as few modes as as possible rather than a continuum is to 
explore the variety of dynamical behavior that can arise 
in a non-dissipative system and its sensitivity to initial 
state. 

In order to better qualify the entanglement dynamics, 
we must specify a way of quantifying entanglement. For 
two qubits one may compute the entanglement of for- 
mation efficiently in terms of the Wootters concurrence, 
and the concurrence itself may be used as an entangle- 
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merit monotone This is the quantity we will use 

most often to quantify the atom-atom entanglement in 
this paper. Alternatively, to quantify entanglement one 
can also use the negativity [32j (which quantifies the de- 
gree to which the state violates the Peres-Horodecki pos- 
itive partial transpose condition). This is useful on larger 
Hilbert spaces, and we use this to quantify entanglement 
when considering aspects of the the entanglement dy- 
namics involving field states. 

In section |TI| we describe a useful way to parameter- 
ize the class of Hamiltonians that can give rise to dif- 
ferent entanglement dynamics due to atomic separation. 



In section [IV| we focus on two special cases in terms of 
separation with the aim of illustrating the variety of dif- 
ferent behaviors that can arise in different initial atomic 
and field states. Specifically, we will show how SD, dy- 
namical entanglement generation, and other phenomena 
differ between these two cases by comparing with anal- 
ogous situations of the two well-studied types of models 
- where each system is a single two-level atom interact- 
ing with a single field mode [TTJ |12l [23] and two two- 
level atoms interacting with a single common field mode 
[24 , 27J. We will show that the entanglement dynamics in 
these models differs significantly to illustrate the effect of 
atomic separation. Finally, in section [V| we examine the 
entanglement dynamics for other atomic separations, us- 
ing analytic results for a few cases and numerical results 
more generally, to examine to what degree the behaviors 
of the special cases examined generalize. To look at the 
entanglement evolution for intermediate values of separa- 
tion analytically, we use the resolvent expansion method 
of [20j restricted to initial states with a single excitation 
among the atoms and field modes. Also, for the most 
general case we can look at the dynamics by numerical 
diagonalization of the Hamiltonian for all the other initial 
states and arbitrary separations. In this way we explore 
the implications of interatomic distance on entanglement 
dynamics in this simplest non-trivial context. 



II. SYSTEM HAMILTONIAN 

Consider a pair of identical two-level atoms coupled 
via the minimal coupling interaction Hamiltonian to a 
collection of electromagnetic (EM) field modes under the 
dipole and rotating-wave approximation. We will assume 
that atoms can be effectively considered to have fixed 
center of mass positions Rj . The interaction Hamiltonian 
of the system in the Schrodinger picture has the form 



where aj and are the raising and lowering operators 
for the two-level system representing the j*^ atom, and 
dl „ and CLq^s are the creation and annihilation operators 
corresponding to the q^^ normal mode with frequency Uq 
and polarization indexed by s in some mode decomposi- 
tion of the EM field with (complex) classical electric field 



mode function /^^ (r). Let us assume that the atomic 
transition in question does not alter the angular momen- 
tum of the atom. Under these assumptions 
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(2) 



where the field is quantized in a region of volume V, dj 
is the complex dipole matrix element for the transition 
in the j*^ atom, and the dipole moment has a direction 

In general, the mode functions /^^ {Rj) can be quite 
complicated, and even for a single atom there can be 
position dependence in the dynamics arising from the 
boundary conditions that the mode functions obey. In 
order to distinguish these effects from the effect of atomic 
separation, we will consider mode functions of the form 



fq V) 



(3) 



in terms of some set of coordinates {u, w), with 5 = 1,2 
indexing the polarization. This form describes traveling- 
wave mode functions that are solutions to an EM bound- 
ary value problem which is translation invariant under 
translations in the coordinate and it appears in a 
number of interesting physical systems, including spher- 
ical resonators [25] and toroidal resonators [26, 38j. We 
will consider situations in which all atoms share the same 
coordinates u = uq and v = vq, differing only in the co- 
ordinate Wj^ so that the variation in behavior due to the 
atomic separation will manifest without additional posi- 
tion dependence arising from the boundary conditions. 
Let us further assume that the atoms are arranged such 
that the atomic transition dipole vectors are aligned with 
one of the mode polarizations, so Cdj ■^q,s = ^s,i- In this 
situation only one polarization is relevant, and all po- 
larization labels will be suppressed. With these further 
assumptions we may express the coupling constants as 



gj^q = -idjfq {uo,vo)e 



2hwqeoV^ 



(4) 



however, it turns out that without loss of generality one 
may study a much smaller set of Hamiltonians. 

Given a particular set of complex phases for the cou- 
pling constants gj^q^ one may make a trivial basis trans- 
formation 
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which simply amounts to redefining the reference for the 
phases of the atoms by cF^e~^^^ and for the fields 

by dq dqe~'^^'i . After this transformation, the interac- 
tion Hamiltonian becomes 

H'j = U.HjUl = hY^Yl djA^-^ + (9-,,) * (6) 



with gj^^ = gj^qC *^^e and we may obtain the so- 
lution to the original problem by using the transformed 
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Hamiltonian: 

\^ (t)) =e-^^^^/^ \^ (0)) = Uy^'^'^^U, IV^ (0)) 

=fjtg-i^W^|^^(0)). (7) 

In terms of the entanglement dynamics, we are not con- 
cerned with the basis transformation as it only acts lo- 
cally on each of the atoms and field modes. If we con- 
sider two identical atoms, with frequencies ujq and dipole 
strengths \dj \ = d coupled to a single EM field mode, then 
the previous paragraph implies that it suffices to consider 
only a Hamiltonian where both couplings are real, and 
the atomic separation does not enter; thus, there can be 
no non-trivial distance dependence. However, if we con- 
sider the two identical atoms coupled to two field modes, 
then without loss of generality we can write the total 
Hamiltonian of the system as 



+ {gi&^di 



+ ujia[di 
-gicr^di 4 



(8) 



where all distance dependence arises from (j) = 
{k2 — ki) {w2 — wi), and "H.C." denotes the Hermitian 
conjugate terms. For a pair of counter-propagating 
modes, we have k2 = —ki = 27r/A and hence the pa- 
rameter (j) = AirAw/X. For simplicity, we will further 
assume that uji = UO2 = ^o? which implies gi = g2 = g. 

Our aim is to get a sense of the variety of different en- 
tanglement dynamics that can result from different sep- 
arations and initial conditions. To that end, we will 
first consider two special cases, which are arguably ex- 
treme cases of the general model, (j) mod 27r = (where 
Aw = nX/2) which we will call the case of symmetric cou- 
pling, and (j) mod 27r = tt (where Aw = (1/4 + n/2) A), 
which we will call the case of antisymmetric coupling (of 
the second atom to the second mode). Next we will dis- 
cuss how the qualitative features of entanglement dynam- 
ics can differ significantly between these two special cases 
of atomic separation. 



III. QUANTIFYING ENTANGLEMENT 

In order to quantitatively study the entanglement dy- 
namics, we must specify a way of quantifying entangle- 
ment. For two qubits one may compute the entangle- 
ment measure known as the entanglement of formation 
efficiently in terms of the Wootters concurrence [31j. In 
order to compute the concurrence, one must choose an 
arbitrary basis and construct the spin-flipped operator 
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where the star denotes the complex conjugation in that 
basis and cry is the Pauli matrix in the same basis. The 
concurrence is 

C (p) = max {0, v/aT - v/A^ - v/A^ - 7^1} , (10) 



where the XjS are the eigenvalues of the matrix pp in 
decreasing order. The entanglement of formation can 
then be computed as 



Sf (C) = h 



1 + Vl - C2 



where 



h (x) = -X log2 {x) - (1 - x) log2 {l-x). 



(11) 
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Alternatively, to quantify entanglement one can also use 
the negativity [32j (which quantifies the degree to which 
the state violates the Peres-Horodecki positive partial 
transpose condition), defined 



M{p) 



(13) 



where p^^ is the partial transpose of the density matrix, 
defined by 



{jMp^^\i,q) = {i,q\p\i,k) 



(14) 



in an arbitrary basis B, and ||-||^ is the trace norm. Either 
of these approaches can be used to quantify entangle- 
ment. The negativity is arguably less precise, because in 
some systems it can be zero even for an entangled state, 
but we will not consider such situation. The advantage 
of the negativity is that it can be efficiently computed 
even in a large Hilbert space. We will most often use the 
concurrence in this article; however, we do quantify the 
entanglement of field states using the negativity. 



IV. SPECIAL CASES 

As a first step toward evaluating the influence of in- 
teratomic separation on entanglement dynamics, we now 
examine two particularly simple cases that one might ex- 
pect to represent two distinct extremes. Since we have 
assumed the two field modes have the same frequency, 
rather than using the original modes Fi and F2 one could 
equally well choose a different mode decomposition of the 
field where the two modes are replaced by linear combi- 
nations TFi and TF2 with annihilation operators 



1 

Ai = — (ai +a2) 
A2 = (ai - 0.2) • 



(15) 
(16) 



These will be two orthogonal standing wave modes. The 
assumed translation invariance of the mode structure in 
the direction of the coordinate w allows us to choose the 
phase of the standing wave modes arbitrarily (subject to 
the constraint of orthogonality), and the mathematical 
definitions we have given here specify that TFi will have 
an anti-node at the position of atom 1 and TF2 will have 
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a node at that location. The couphng to atom 2 then de- 
pends on the separation. With this mode transformation 
we can now turn to the special cases of interest. 



When = we have the case of two modes with sym- 
metrical coupling (SC) to the second atom. The interac- 
tion Hamiltonian can then be written 



Hi = V2hg aj^Ai + + cr^ M + cr^Al (17) 



in terms of these transformed modes, so that the atoms 
only couple to Ai and not 74.2. This reflects the fact 
that the atoms are separated by Ak; = nA/2, so the sec- 
ond atom also lies at a node of TF2 and anti-node of 
TFi. Thus, the Hamiltonian is equivalent to the model 
of a single mode symmetrically coupled to two atoms, a 
Tavis-Cummings (TC) model, and all work done on that 
model can be applied here. The evolution of the reduced 
density matrix of the atom pA (t) in the SC model is 
be the same as in the TC model with a properly trans- 
formed initial state. Namely, if the total system (atoms 
and modes) is in an initial state described by the density 
matrix xsc (0), then the appropriate initial density ma- 
trix for the TC model is obtained by making the mode 
transformation of Eqs. 



(15 



and (16) and tracing out 
the second transformed mode TF2. As this mapping of 
initial states involves a partial trace, it is a many- to-one 
mapping from the SC problem to the TC problem (and 
this mapping does not preserve purity). 



In order to solve the dynamics in the TC model, and 
by extension the SC model, we simply compute the 
time evolution operator expressed in the atomic basis 
{|ee) , \eg) , \ge) , 1^'^)} directly by exponentiation (as in, 
e.g. 1221): 



Hi=V2g 
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(18) 



where 



Si=A 



C2 = 





COS (^y/iAgtj 
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A 
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(19) 



C3 = \ (cos (Vilgt) +1) Ci = ^ (cos (^Vilgtj - 1) 



A,^A, 



Ce = l-A\^A,^ 
A = AiA\ + A\Ai 



cos (\AAqt] 



cos 



A 



Ai 



We observe that there is generally oscillatory behavior for 
the state dynamics coming from the periodic time depen- 
dence and the frequency of these oscillations depends on 
the number of excitations in the initial field state. Also, 
a Fock state has a well defined frequency. 

When ^ = TT we have a two-mode model with asym- 
metric coupling (AC) to the second atom. In terms of the 
transformed modes of Eqs. ( 15 ) and ( 16 ) the Hamiltonian 
can be written as 



Hi = V2hg a^Ai + + a^A2 + A^ 



(20) 



In this case, rather than both atoms coupling to one 
mode we see that Atom 1 couples only to transformed 
mode TFi while Atom 2 couples only to TF2. Because 
in this case Aw = (l/4 + n/2)A, atom 2 now lies at a 
node of TFi and anti-node of TF2. This situation is 
equivalent to a model comprised of two subsystems that 
are totally isolated from one another, each composed of 
a single atom coupled to a single mode. We will call this 
the double Jaynes-Cummings (DJC) model. This sort of 
model with isolated subsystems is common to the study 
of entanglement sudden death and the DJC model 

specifically has been studied pT| \T2\ |23] , so the mapping 
allows those investigations to bear on this model. 

As in the previous case, the evolution of pA {t) for the 
AC model should be the same as given by the DJC model 
with the proper mapping of initial states. This implies 
that there can, for example, be no dynamical genera- 
tion of atomic entanglement in the AC model unless the 
DJC initial field state obtained by the mapping is entan- 
gled. In the DJC model we can write the unitary time- 
evolution operators for the two separate non-interacting 
atom- field subsystems as /7i and /72, and then the total 
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time evolution operator is U = /7i (g) We again com- 
pute the two subsystem unitary time evolution operators 
by direct exponentiation to obtain 



U2=e- 
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in applying it only to instances where the entanglement 
goes to zero for some finite (non-zero) duration. In the 
case where entanglement goes to zero only for an instant 
during the time evolution we refer to it as death for an 
instant (DI). If there is a non-zero entanglement at all 
times once it is generated in the system then we label 
it as being "always living" (AL). In the context of that 
general summary, we now turn to commenting on some 
specific features. 



Entanglement generation in the AC model 
Transfer of correlations 



with 



Ca = cos{ J 2A,Algt) Q2 = cos{ J 21} Aigt) (23) 



_sm{^2AAlgt) _ sm{\/2AlA,gt) 

^il — / Oi2 — A - — 



In order to illustrate the variety of behavior that 
can arise among these models, we examine the en- 
tanglement dynamics of a selection of initial states 
in which the atoms are uncorrelated with the fields. 
In each case we select the atomic state from the set 
of pure states {|^^) , |ee) , |e^) , |^) , |^)}, where |$) = 
(|ee) + 1^^)) / A 1^) ^ {\eg) + \ge)) /x/2. The ini- 
tial field state will be either a product of Fock states 
|nAr,mAr), a product of Glauber coherent states |ac,/3c), 
a product of squeezed vacuum states \isq-, ~isq) 1 a two- 
mode squeezed vacuum state (TMSS) |^, 0, 0) , a thermal 
state (g) p^^^ (with both modes having equal tempera- 
ture), the pure state \rinm)i oi" the mixed state pum- The 
two-mode squeezed state is defined as the state result- 
ing from the action of the two-mode squeezing operator 
S{£) = e(^*^i^2-^^^^2) on vacuum. The state \rjnm) is 
the result of mapping the state |nAr,mAr) in the origi- 
nal modes of the AC problem to the transformed modes 
equivalent to the DJC problem, with 



^ n m 

(24) 

y^{k^l)\{-iy \m^n-k-l)\k^l), 

and pnm = T^^TF2[\Vnm) {Vnm\] IS the statc in the 
TC model that gives equivalent evolution to the state 
|^Ar,^Ar) in the SC model. To find the atom-atom den- 
sity matrix we trace out the field modes from the time 
evolved state of the closed atom-field system. 

We summarize our findings for the entanglement be- 
havior given the various initial states considered in the 
four models in Table |lj listing the equivalent SC-TC cases 
and AC-DJC cases. When discussing entanglement sud- 
den death, we adopt the usage of the term as in [i2\ 



By reducing the SC and AC models to the TC and DJC 
models, respectively, we see that the change in atomic 
separation leads to very different entanglement dynam- 
ics. In the case of the AC model, we have a mapping 
to the dynamics of the DJC model, where each atom in- 
teracts only with a separate field mode. So in the DJC 
model, the dynamics cannot increase entanglement be- 
tween the two atom-field subsystems; therefore, if the 
atomic state is not entangled initially then it will remain 
separable, unless there is an initial entanglement between 
the field modes that can be transferred to the atoms by 
the dynamics. Knowing this, we can see that any initial 
field state for the AC model that maps to a separable 
DJC field state fails to generate entanglement. The na- 
ture of the mapping means that even some entangled 
field states will fail to generate atomic entanglement in 
the AC model, while some separable states will map to 
an entangled DJC state and will generate entanglement. 

For example, when considering the AC model with an 
initial field state that is a product of squeezed states 
l^sg, —^sq) we find that entanglement is dynamically gen- 
erated as shown in Fig. [T] One can understand this from 
the fact that the state —£,sq) maps to the DJC model 
with an initial TMSS, so the generation of atomic entan- 
glement occurs simply because the dynamics transfers 
the entanglement between the field modes to the atoms. 
Further, we find that there is an optimal squeezing pa- 
rameter value for entanglement generation when the field 
state is close to being a maximally entangled qubit state 
(ao(0 |00)+<^i(0 l^))' with ai(^) and ao(0 being com- 
parable. On increasing the squeezing parameter further 
there are contributions from higher Fock states which de- 
crease the transfer of entanglement. 

Transfer of field- field (TF1-TF2) correlations can be 
explicitly seen in a small squeezing approximation when 
a two mode squeezed state interacts with an initially 
separable atomic state |ee). Since the probability of 
higher photon numbers is small, by restricting to the 4- 
dimensional subspace of lowest energy states in the Fock 
basis for TF1-TF2 we obtain the negativities for the two 
subsystems as. 
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Initial field State 


Atomic State 


SC 




A 1 \ 

A. \ee) 


B. \eg) 


C. \gg) 


D. 1$) 


E. 1^) 


1. \nN,mN) 
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No 


Yes, DI 
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DI/SD 


SD 


SD 






No 
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Yes, SD 


SD 


SD 


3. , 


Oth 




No 


Yes, DI 


Yes, SD 


AL/SD 


SD 


4. 


^(«= + /3c),^(ac-/3c)) 




Yes, AL/SD 


Yes, SD 


Yes, AL/SD 


AL/SD 


AL/SD 


5. 






No 


Yes, DI 


Yes, SD 


AL/SD 


SD 


6. 1^,0,0) 


16.) 


Yes, AL/SD 


Yes,SD 


Yes, AL/SD 


AL 


SD 




Initial field State 


Atomic State 


AC 


DJC 


A. ee) 


B. \eg) 


C. \gg) 


D. 1$) 


E. 1^) 


1. \nN,mN) 


1 '^nm ) 


Ye^^ SD 


Ye^^ SD 


Yei 
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SD/DI 


SD 


SD/AL 




njv,miv) 


No 


No 


No 


SD 


SD 


3. Pth 




No 


No 


No 


SD 


SD 


4. |ac,/3c) 




No 


No 


No 


SD 


SD 




e,o,o) 


Yes, SD 


Yes, SD 


Yes, SD 


SD 


SD 


6. 1^,0,0) 


— 6g) 


No 


No 


No 


SD 


SD 



^ No entanglement for = ttt^n 

Table I: Entanglement dynamics for two modes symmetrical coupling (SC) with </) = and anti-symmetrical 
coupling with = tt. Columns A-C list whether there is an entanglement generation in an initially separable atomic 
state (yes or no). The dynamical phenomena observed in columns A-E are listed as entanglement sudden death 
(SD), entanglement dies for only an instant (DI), entanglement remains non-zero at all times and so is "always 
living" (AL). A '/' denotes that both kinds of dynamics are present depending on the particular initial state chosen 

from within the class indicated. 



where si = sm{\^gt), ci = cos{^/2gt) and C2 = cos{2gt). 
Fig. [2] shows how the dynamics of the system transfer the 
entanglement between field-field subsystem and atom- 
atom repeatedly. 

B. Entanglement generation in the TC model and 
initial states with identical entanglement dynamics 

Turning to the SC model, the correspondence drawn 
to the TC model leads to the essential feature that the 
map for the initial field states from SC to TC is many 
to one; the set of initial states that have the same re- 
duced density matrix for the first transformed mode TFi 
have identical entanglement dynamics in terms of atom- 
atom entanglement. As a counterintuitive example of 
this feature we note that a squeezed state of the form 
\£,sq^ ~^sq) and a thermal field in the SC model give the 
same entanglement dynamics provided the average num- 
ber of photons in the thermal field corresponds to that 
in the squeezed state {nth = sinh^ l^sq])- 

We can also turn to previous work on the TC model 
to expose interesting features of the SC model. As had 
been shown in [24j, in the TC model an atomic state 
|ee) interacting with any Fock state |n) in the field mode 
never gets entangled. This can be generalized to any ini- 
tial field state with a density matrix that is diagonal in 
the Fock basis as follows: If one defines the time evolved 



state It/jn) = U \ ee) |n), then the atom- atom density ma- 
trix given as p'^^^ = Tip {^pn\] remains separable. 
Extending to a general field density matrix diagonal in 
the Fock basis, the time evolved atom-atom density ma- 
trix given as pi2 = Tip [}Zn -^^P^^"^] is clearly a convex 
sum of separable density matrices, and hence there is no 
entanglement generation. 

Looking at the initial state \eg) \n) in the TC case, 
we observe no SD in Fig. |3a| This can be explained 
by considering the state as a superposition \eg) = 
I {\eg) + l^e)) -\- ^ {\eg) — |^e)), where due to the sym- 
metry of the coupling constants the maximally entangled 
dark state {\eg) — \ge)) />/2 does not interact with the 
field. It is only momentarily during the evolution that 
the state of the system returns to being the original sep- 
arable superposition \eg). As a result we always observe 
some entanglement between the two atoms for an initial 
field density matrix diagonal in the Fock basis. The dy- 
namics of the \gg) \n) state in the TC model, shown in 
Fig. |3b[ exhibit SD in general except for the special case 
of n = 1 where because of symmetry reasons the state os- 
cillates between the states \gg) |1) and |^) |0), going from 
being separable to maximally entangled. Hence, any den- 
sity matrix of the field modes diagonal in the Fock basis 
with a high component of |1) (1| would generate more 
entanglement in general. 

Apart from the Fock state in the TC model and a ther- 
mal field, another example of having a initial field density 
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(a) |ee> (b) \eg) 

Figure 1: Entanglement dynamics for a two-mode squeezed state (|^,0,0)) in DJC or product of single mode 
squeezed states —^5,0)) in AC interacting with different initial atomic states. Entanglement generation from 
transfer of correlations from the field, optimal squeezing for maximum transfer of correlations can be seen. 



matrix 




Figure 2: Small squeezing entanglement dynamics for a 
two-mode squeezed state in the DJC model interacting 

with an initially separable atomic state |ee), 
entanglement is generated in the atomic subsystem via 
transfer of field-field correlations. The red dashed line is 
for the field-field entanglement, blue solid line is the low 
temperature approximation of the atom-atom 
entanglement and black dashed is for the exact 
atom-atom entanglement 



matrix diagonal in the Fock basis is to have a Fock state 
|n, m) in the SC which corresponds to the TC density 



n m 



Pnm — ^ ^ 



f^mnki \m -\- n — k — I) {m -\- n — k — l\ 



2^+^n!m! 



(26) 



where i^mnki = ''Ck''Cp'^CrCqSk+i,p+q{m ^ n - k - 
l)\{k + /)!(— 1)^ For this state again we see no entan- 
glement generation for |ee), DI for \eg) and maximal en- 
tanglement in \gg) |10), as shown in Figs. 3e and[3f[ 

Another point we observe is that for an initial state 
|ee) |n,n) there is no entanglement generation, which is 
a common feature between SC and AC. A Fock state 
|n,m) in AC transforms into an entangled state in the 
DJC model, so we expect to have entanglement gen- 
eration in the system for an initially separable atomic 
state. As an exception we find that for n = m, if there 
is no atom-atom entanglement to begin with then the 
atoms remain separable evolving to a diagonal density 
matrix. This is counterintuitive in the sense that there 
is no entanglement transfer to the atomic subsystem 
from the entangled field modes. This feature can be 
explained by considering the initial field state \irinn) = 
2h. ^k=o "Cfc(-l) V2^K2n - 2k)\ \2n - 2k, 2k). If the 
atoms are initially in the state \ee) then time evolution 
will lead to an entanglement of atomic and field states 
such that detecting whether the number of photons in 
each mode is even or odd tells us the state of the two 
atoms. On tracing out the field this gives us a diagonal 
atom-atom density matrix with no atom-atom entangle- 
ment. The same is true for the initial atomic states \eg) 
and \gg). 
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(e) \eg) with |niv,^Ar} in SC or pnm in TC (f) \gg) with |niv,'^Ar} in SC or pnm in TC 

Figure 3: Entanglement dynamics in presence of initial field states diagonal in the Fock basis of the TC - no 
entanglement is seen for initial atomic state |ee), only DI is observed for the atomic state \eg) and maximum 

entanglement for n ~ 1 for |^^) 
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C. Selected contrasts in entanglement dynamics 

While Table [T| encapsulates the qualitative features of 
entanglement dynamics for many states, it may be use- 
ful to point to a few interesting examples. Note that 
many familiar initial field states that are separable in 
the DJC model fail to dynamically generate entangle- 
ment in the AC model. In contrast, entanglement gener- 
ation is a common feature in the SC model for the same 
selection of field states. For example, starting with a two- 
mode squeezed state in the field modes does not generate 
any entanglement in AC, because the corresponding field 
state in DJC possesses no correlations amongst TFi and 
TF2. On the other hand, a two- mode squeezed state in 
the symmetrically coupled case generates entanglement 
in the initially separable atomic states as can be seen 
in Fig. [4j Furthermore, for the separable initial atomic 
states \ee) and \gg)^ if the field is sufficiently squeezed, 
entanglement is dynamically generated and once gener- 
ated sustains forever, protected from sudden death. This 
behavior is shown in Figs. 4a and 4c If the atoms are 



initially in the entangled atomic state |<l>), then Fig. 4d 
shows that the entanglement stays protected as the state 
becomes almost a dark state for large squeezing param- 
eter values. One can also see entanglement protection 
for a thermal field in the SC model so long as the ther- 
mal average photon number is below a threshold value 
ficrit ^ 0-43 (Fig.[5| with the state being AL (above this 
critical temperature |^) experiences SD as well). 

In considering a selection of familiar states for these 
two special cases of symmetric and anti-symmetric cou- 
pling, we find that the same sort of initial states can have 
very different entanglement dynamics when the spac- 
ing of the atoms is altered. For the selection of ini- 
tial states considered, the symmetric coupling, where 
/^w = nA/2, is better for generating entanglement be- 
tween the atoms, and for protecting entanglement from 
sudden death. However, these remain only special cases 
that reveal some facets of the differences in entanglement 
dynamics arising from a change in distance. It remains 
to be seen what happens for a more general selection of 
distances. 



V. GENERAL SEPARATION 

To extend the results found in the previous section 
to arbitrary atomic separations we use two different ap- 
proaches to solve the problem — the resolvent approach 
as described in [20] and numerical diagonalization of the 
Hamiltonian. The first approach gives an analytic so- 
lution but only yields sufficiently simple expressions to 
be of use for initial states that have at most one exci- 
tation among the atoms and field modes, meaning su- 
perpositions and mixtures of states belonging to the set 
= {|e^00) , l^eOO) , \ggl^) , \gg^l)}. Using the re- 



solvent expansion of the Hamiltonian 



iHt 



(27) 



and expanding 



{E - H)-^ = {E- Ho)-^ + (^ - Ho)-^Hi{E - Ho)'^ 
+ (E - Ho)-^Hi{E - Ho)-^Hi{E - Ho)'^ + ... (28) 

we can find an exact expression for the matrix ele- 
ments Tij = \ {E — H)~^ by resumming the se- 



ries expansion in ( |28| ). It can be seen that the action 
of the operator A = {{E — Hq)~^Hi)'^ on the subsets 
= {|e^OO),|^eOO)} and = { 1^5^10) , |^^01)} 

gives a state in the same subspace. Say, A{jpn^ |<^i) + 



<ln^ l«2» = Pnii + Ci l«2) and A{p\t' 



(1) 



.(2) 



<ln^ 1/52)) = P'nii + q'nli \p2) , SUch that 

/ (1.2)\ / (1,2)\ 

then, 



.(2) 



J2) 



(29) 



2ff^ / 1 e'^/2 cos((/,/2) 

{E-ujqY U-^*/2cos(0/2) 1 



(30) 

and A^'^^ = A^^^ . Now consider the matrix element 
^11 = (e^00| {E — H)~^ \egOO)^ by using the expansion 
in ( 28 ) this can be summed up to obtain 



Til 



1 



1 



1 

4 V^-^o -2^sin(0/4) ' - cjq + 2^ sin(^/4) 
1 1 



E -ujQ-2g cos((/)/4) E -ujQ^2g cos((/)/4) 

: T22 = T33 = T44 



(31) 



Similarly, 



-12 



4 - cjo - 2^cos((/)/4) £; - cjo + 2^cos((/)/4) 

1 1 



E -uJo-2g sin(^/4) E - ujq ^ 2g sin(0/4) ^ 
: T*i = T3*4 = T43 (32) 

g-i0/4 /I I 



-13 



4 - ^0 - 2^cos(^/4) - cjo + 2^cos(^/4) 

i i 



E -uJo-2g sin((/)/4) E -uoQ^2g sin((/)/4) 



^3*1 



23 



^32 



^246 



-i(f) 



-icf) 



(33) 
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1 2 3 4 5 6 

gt 

(c) \gg) - SD disappears on increasing the squeezing 




1 2 3 4 5 6 

gt 

(d) \^) - State moves towards being maximally entangled at all 
times as squeezing is increased 



Figure 4: Entanglement dynamics for a single mode squeezed state (|^s)) in TC or two- mode squeezed state 
(1^,0,0)) in SC interacting with different initial atomic states 



From those matrix elements we obtain the matrix ele- 



ments of the propagator as in (27): 



U21 



U:. 



31 



(cosi^it + cosi^2^) = U22 = U33 = Um 



{cos Kit — cosi^2^) = U4S 



{cos Kit — cosK2t) = U'^ 



34 



{—i sin Kit + sin K2t) = [/41 



C/32 = C^24e-'' 



{—ismKit — sini^2^) = ^14 



^23 ^426' 



i0 



(34) 
(35) 

(36) 

(37) 
(38) 



where Ki = 2g sin | and K2 = 2g cos |. 

It can be seen as a general observation that for any 
initial state that lies in the single excitation subspace 
the time evolved atom-atom density matrix has the form 

P = Pu \eg) {eg\ + P22 \ge) {ge\ + pi2 \eg) {ge\ 

+ P12 l^e) {eg\ + pss \gg) {gg\ (39) 

[22J. For this general form of the density matrix the 
atom-atom concurrence simplifies to C = 2|pi2|, meaning 
for a state in the single excitation subspace the entangle- 
ment of the two atoms has a dynamics that is either AL 
or DI (since pi2 is analytic and can never vanish over a 
finite non-zero length of time). 
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Figure 5: Entanglement dynamics for a single mode 
thermal field interacting with an initially entangled 
atomic state |^) in the SC model. SD occurs after a 
certain threshold temperature. 



Using the propagator for the one-excitation subspace, 
we can examine the analytic expressions for entanglement 
dynamics of a few initial states. We begin with the state 
\egOO). The time evolved reduced density matrix of the 
two atoms after tracing out the field modes is obtained 
as 

Pi2 = \Unf \eg) {eg\ + \U2if \ge) {ge\ + UnU^^ \eg) {ge\ 
+ U^,U2i \ge) {eg\ + (IC/31I' + IC/41I') \gg) {gg\ (40) 

This gives a concurrence for the two qubits as 

C{(t),t) = i |cos(2i^it) - cos(2i^2t)| (41) 

It can be seen from the above expression that, unless (j) = 
TT, the atoms are entangled except at a discrete, periodic 
set of instants given by td = m7r/{Ki — mel where 
the state is separable (i.e., "dead"). As the separation 
parameter (j) approaches tt, Ki approaches K2 and the 
death time approaches infinity. Thus for antisymmetric 
coupling, there is no entanglement generated at all, as the 
situation corresponds to having the two atoms interacting 
with two separate transformed modes. Thus, this special 
case is uncharacteristic of the qualitative features of the 
entanglement dynamics for general separations, and for 
any other separation the atoms will become entangled, if 
perhaps on a long timescale. 

For the initial state \gglO) the atom- atom entangle- 
ment is 

C{cP,t) = i {\sm\Kit)\ + \sm\K2t)\) . (42) 

In this situation it can be seen that the entanglement 
remains AL except for separations where tan((/)/4) is 



rational, where the entanglement dynamics is DL This 
latter case includes symmetric and antisymmeric cou- 
plings, where the concurrence oscillates as C(0,t) = 
l/2|sin^(2^t)| and C(7r,t) = |sin^(V^^t)| respectively. 
Here we see that the special cases describing states with 
high symmetry are not at all representative of the generic 
qualitative features for arbitrary separation. For the ini- 
tially entangled atomic state ^(|e^) + \9^)) I^O)? ^^e 
concurrence is given as C{(j)^t) = cos^(i^it) cos^((/)/4) + 
cos^(i^2^) sin^(0/4), again we can see that entanglement 
is AL unless tan(^/4) is a rational number. 

Restricting to the single excitation subspace, one might 
also seek to find the approximate concurrence for thermal 
and coherent states with low average photon numbers 
(keeping only one excitation in the field modes) interact- 
ing with \gg); however, one can find that the resulting 
values are of the same order as two-excitation terms that 
have been neglected, so this approach is not useful in that 
case. 

To look at the entanglement dynamics for other ini- 
tial states we numerically diagonalize the Hamiltonian, 
restricting to a small photon number approximation. We 
find that SD is generally observed for most distances 
with the initial states we have discussed in the preced- 
ing section. In many cases, however, the special cases 
of symmetric and antisymmetric coupling show qualita- 
tively different dynamical behavior. For example, if we 
look at the entanglement dynamics for the maximally en- 
tangled state \^) interacting with a thermal state in the 
SC model, we had seen that the entanglement was pro- 
tected upto some critical temperature, while there was 
SD for the AC model. Upon extending this to arbitrary 
separation values we find that for all separations where 
7^ the entanglement always shows SD for all temper- 
atures as shown in Fig. [9j Similarly, for |^) interacting 
with a two-mode squeezed state one has SD except for 
symmetric coupling or (/) = 0. As we had discussed in the 
previous section for the case of an initial EM field state 
diagonal in the Fock basis in the TC model, if both the 
atoms are excited there is no entanglement generation 
and for the atoms in \eg) the entanglement dies only for 
an instant. Considering the Fock state |nAr,rZiv) in the 
SC model translates to the state pnn in the TC model 
which is diagonal in the Fock basis, hence we see no en- 
tanglement generation for \ee) and DI for \eg). On the 
other hand, we had also shown that for such a state there 
is no entanglement generation in the AC situation. We 
use numerical approach to consider the entanglement be- 
havior for intermediate separations. It can be seen from 



Fig{TO]that the intermediate distances have sudden death 
for the case of \egnn), even if one deviates slightly from 
the symmetric coupling (j) = 0. The state \eenn) shows 
entanglement generation for asymmetric couplings. 

While examining the special cases with high symmetry 
in the setup it is useful to show the over-riding effect of 
symmetry on the entanglement dynamics of this system, 
computing the results for arbitrary distances suggests 
that these special cases are not representative in terms of 
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the general qualitative behavior of entanglement. As we 
have discussed above, in many cases the dynamics falls in 
one qualitative class almost everywhere on the interval of 
asymmetrically coupled cases (0 / or tt) but in a dif- 
ferent class on one or both boundaries. Moreover, of the 
three qualitative classes we have discussed here, all are 
possible at every distance, given the correct initial state. 
An analysis of linear oscillators coupled to a common 
vacuum field (consisting of a continuum of modes) by 
Lin and Hu found that there exists an entanglement dis- 
tance, so that oscillators with a separation less than that 
distance will be asymptotically entangled, while initially 
entangled oscillators spaced further will undergo SD [35j . 
In contrast, an analysis of two two- level atoms interacting 
with a common EM thermal field (with a continuum of 
modes) found [36j that all initial states will eventually un- 
dergo SD for any non-zero separation. In this two-mode 



model we find no evidence of an "entanglement distance" 
which determines a change in the qualitative features of 
entanglement dynamics. For each of initial states consid- 
ered, the qualitative class into which the entanglement 
dynamics falls is the same for almost all distances, with 
departures only on a set of measure zero. For a thermal 
field, which is the environment considered in [35',"36] (the 
former considers specifically the vacuum), SD occurs at 
almost all distances. This bears some resemblance to the 
result of [36], though the non-dissipative nature of this 
model ensures revival and, in this regard, is qualitatively 
quite different. 
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(a) (b) 
Figure 8: Entanglement dynamics for -^(\eg) + \ge)) |00), AL/DI behavior for all values of separations 





(a) Initial state |$} interacting with a thermal state with 
average number of photons nth = 0-2 



(b) 1$) interacting with a two-mode squeezed state with 
6a = 0.1 



Figure 9: Entanglement dynamics for |<l>), AL behavior for symmetric coupling, SD for all the other separations 



VI. SUMMARY AND DISCUSSION 

We have analyzed the entanglement dynamics in a 
model consisting of two two-level atoms and two electro- 
magnetic field modes for a variety of familiar field states 
and classified the various cases in terms of different dy- 
namical behavior. One aim of this analysis is to get a 
sense of the variety of entanglement dynamics that can 
arise from different atomic separations in a case where 
two atoms interact with a shared EM field. It is useful to 
examine this question in a very simple model that can be 
solved with a minimum of approximations (using which 
can lead to unphysical effects) and understood in detail. 
We have argued that there is no non-trivial distance de- 
pendence in a single mode model, and, therefore, a two- 
mode model represents the simplest case for the study of 



distance dependence in the entanglement dynamics. On 
the other hand, including a continuum of modes leads 
to dissipation and misses the variety of dynamical be- 
havior that can be seen with a few number of modes. 
In particular, the effect of initial state is largely washed 
out. We have first studied two arguably extreme cases 
out of the class of two- mode Hamiltonians that can arise, 
one in which the two modes are symmetrically coupled 
(SC) to the second atom, and one where the two modes 
have asymmetric coupling (AC). A useful insight in un- 
derstanding these models is that the atomic dynamics 
in the SC model correspond exactly to the dynamics for 
a Tavis-Cummings model with a single field mode sym- 
metrically coupled (TC) to both atoms with a suitable 
mapping of the field state, while the atomic dynamics 
for the AC model correspond exactly to the dynamics 
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(a) \eg22) has DI for symmetric coupling {(j) = 0) and no (b) |eell} shows no entanglement generation for the 

entanglement for antisymmetric coupling = tt), SD for all symmetric and anti-symmetric coupling, SD otherwise 
the intermediate separations 



Figure 10: Entanglement dynamics for Fock states with general separations 

for a double Jaynes-Cummings (DJC) model, made up thermal generation of entanglement in [27J. We extend 
of two isolated subsystems each with one atom coupled the conclusion of these results more generally to the set 
to one mode, under a proper transformation of the field of initial field states that are diagonal in the Fock basis, 
state. These mappings help one understand the signifi- We find that for a field state diagonal in the Fock basis 
cant differences in behavior in the two seemingly similar there is no entanglement generation in |ee), DI for \eg) 
models, giving a window into how symmetry introduced and for \gg) an initial field state with a high component 
by specific atomic separation can affect entanglement dy- of |1) (1| would generate more entanglement in general, 
namics. Another significant implication of the mapping Other than that we study the case of an initially entan- 
between the SC and TC models comes from the fact that gled atomic state interacting with a thermal field and 
the mapping of initial field states between those models find that the entanglement is preserved for low tempera- 
is many-to-one, this shows us entire classes of field states tures and the incidence of sudden death occurs only after 
for the SC model that will give exactly the same atomic a certain critical temperature. The DJC model was pre- 
dynamics. viously studied from the viewpoint of transfer of entan- 

In examining the dynamical generation of entangle- glement phenomenon in [HI [T2J and more generally by 

ment from an initially separable atomic state, we find [23] for the case of an entangled atomic state interacting 

quite a marked contrast between the SC and AC models. with a vacuum field and transferring correlations via the 

While entanglement generation is a relatively common atom-field entanglement. We look at the phenomenon 

feature, present for a variety of field states, in the SC for the reversed situation, where one starts with a max- 

model, it is comparatively much more rare in the AC imally correlated field state (TMSS) and looks at the 

model. This difference is not so surprising however, if entanglement generated in the atomic subsystem. Since 

one views it in terms of the mappings we have intro- the field state does not remain gaussian as the system 

duced to the other models. When considering that the evolves, we resort to a low temperature approximation 

AC maps to the DJC model, where the two subsystems in order to quantify the entanglement. We also see that 

are isolated, one would expect entanglement generation the transfer is maximum for an optimal squeezing pa- 

to be relatively rare; it can only exist in cases where the rameter value such that the field state approximates a 

field state is mapped to an entangled state whose entan- maximally entangled qubit state. As a potentially use- 

glement can then be transferred to the atoms. With the ful phenomenon we observe cases where we find that the 

SC, by contrast, we have a mapping to the TC model entanglement always stays alive once generated or stays 

where a single shared field mode can readily introduce protected if already present. In particular, for the case 

entanglement between the two atoms. of a squeezed state interacting with an initially entan- 

These models (TC and DJC) have previously been gled atomic state |<l>) we see that the qubits end up in 

studied for some of the cases that we have looked at as an almost dark state and stay almost maximally entan- 

well. For example, Tessier et al in [24J had also analyzed gled for all times. As a natural question one may ask 

the the initial state \een) and coherent state in the TC what would be the extension of these two extreme cases 

model. Also, interaction of atomic states |ee), \eg) and to intermediate separations and how does the dynamics 

\gg) with a thermal field was studied in the context of on one extreme change into a different behavior on the 
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other. We look at the general separations using the resol- 
vent method for states in the single excitation subspace 
to find the concurrence as a function of separation and 
time. We see that these states can only exhibit AL/DI 
behavior based on the general form of the atom-atom 
density matrix as noted in [22j, in particular it can be 
seen explicitly how the behavior changes for the sym- 
metric and anti-symmetric couplings. We can look at all 
the other cases using numerical analysis and we find that 
for all the cases considered there is no critical separation 
behavior as seen in the continuum cases [35l [36]. For ex- 
ample, situations where entanglement stays protected in 
the symmetric coupling case (thermal state and squeezed 
state) shows SD as soon as there is any asymmetry in the 
coupling. 

We have chosen in this analysis to try to isolate the 
effect of atomic separation on entanglement dynamics 
from the position-dependent effects arising from bound- 
ary conditions. For this purpose, we have assumed in 
Sec. [IT] an atom-field coupling that depends on the co- 
ordinate separating the atoms only by a phase factor. 
For clarity we have supposed our two field modes are 
traveling- wave modes in free space, however a more ex- 
perimentally relevant situation would be a toroidal res- 
onator, where there is rotational invariance in the az- 
imuthal direction that satisfies our requirements. In this 
case the two modes of interest would be two resonant, 
counter-propagating whispering-gallery modes. Because 
strong coupling between an atom and a whispering- 
gallery modes of a microtoroidal resonator has observed 
experimentally [34[ [38 ] , there is the possibility of exper- 
imentally probing quite directly the model we have con- 
sidered. However, in the experimental system there will 
be dissipative dynamics arising from emission into other 
modes outside the resonator, so a detailed comparison 
would require either including the dissipative effects in 
the theoretical model or a restriction to the case of suffi- 
ciently strong coupling to the resonator modes and early 



times that the dissipation could be neglected. 

Our analysis of special cases of the entanglement dy- 
namics arising in two atoms interacting with two modes 
suggests a wide variety of different behaviors can arise, 
with qualitative features of the dynamics changing en- 
tirely between the two special cases considered, even with 
the same initial field state. This suggests that under- 
standing the distance dependence of entanglement dy- 
namics for multiple atoms interacting with a common 
field will be quite important for predicting even the qual- 
itative features that may arise. Furthermore, if one has 
the practical goal of dynamically generating entangle- 
ment or protecting entanglement once generated, the spe- 
cial cases we have considered suggest that the ability to 
achieve these goals will be greatly impacted by the sep- 
aration of the atoms. Having gotten some sense of the 
variation in behavior that can arise, we quantify the en- 
tanglement dynamics over a range of atomic separations 
rather than just special or extreme cases that highlight 
the importance of symmetry factors in entanglement dy- 
namics. 

On this point it would be interesting to compare 
the exact results of the present model with the general 
statements on classifying the entanglement behavior [37J 
based on geometrical and topological properties of quan- 
tum state spaces [4j. 
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